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Point defects in self-assembled crystals, such as vacancies and interstitials, attract each other and 
form stable clusters. This leads to a phase separation between perfect crystalline structures and 
defect conglomerates at low temperatures. We propose a method that allows one to tune the effective 
interactions between point defects from attractive to repulsive by means of external periodic fields. 

In the quantum regime, this allows one to engineer strongly-correlated many-body phases. We 
exemplify the microscopic mechanism by considering dipolar quantum gases of ground state polar 
molecules and weakly bound molecules of strongly magnetic atoms trapped in a weak optical lattice 
in a two-dimensional configuration. By tuning the lattice depth, defect interactions turn repulsive, 
which allows us to deterministically design a novel supersolid phase in the continuum limit. 
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I. INTRODUCTION 

Defects are crucial for the determination of macro¬ 
scopic mechanical, optical, and electronic properties of 
solids One key aspect is the mutual effective in¬ 

teractions between point defects such as vacancies or in¬ 
terstitials, corresponding to the lack or excess of crys¬ 
tal particles, respectively. In self-assembled classical 
crystals the effective interactions between point defects 
is attractive for all combinations of defects (vacancy- 
vacancy, interstitial-interstitial, interstitial-vacancy) in a 
wide range 0H7]. This interaction leads to the formation 
of string-like defect clusters j8]. The mechanism behind 
the attraction of defects is a result of non-linear effects 
in the displacement fields of multiple defects and cannot 
be described as a simple pair interaction within elastic¬ 
ity theory. In the quantum regime, the interaction be¬ 
tween vacancies is even less well understood, and it is 
an open question whether the classical results may be 
directly used to infer many-body properties in the quan¬ 
tum regime. There, the interaction between vacancies is a 
crucial part of the theory for the supersolid phase ehh]. 
The supersolid is conjectured to be a result of delocal¬ 
ized vacancies. However, the precise role of the defect 
dynamics to establish both superfluid and crystalline or¬ 
ders has been the object of intense investigations in the 
last decades p~5I428] . 

Here, we show that effective interactions between point 
defects in a self-assembled crystal can be systematically 
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tuned from attractive to repulsive by means of external 
periodic superlattices. The lattice spacing of the opti¬ 
cal lattice is chosen to be identical or double the lat¬ 
tice spacing of the self assembled crystal. In this setup, 
the interaction between point defects is the result of an 
interplay between displacement induced attraction and 
energy induced repulsion. By properly choosing the pa¬ 
rameters of the superlattice, the relative strength of the 
displacement-induced energy and entropy parts can be 
tuned, thus making the interaction an experimentally ac¬ 
cessible variable. 

Based on this microscopic picture of tunable defects, 
we investigate the many-body quantum regime. As an 
example, using exact quantum Monte-Carlo methods, we 
study the phases of a dipolar bosonic gas in the pres¬ 
ence of a weak triangular optical lattice. Unlike previous 
works that focused on the tight-binding limit valid for 
deep lattice potentials nanaBniEii, here we focus on 
the continuous space limit where the band structure of 
the lattice is not formed or barely formed. We inves¬ 
tigate the phase diagram by varying the lattice depth 
and the strength of the dipole-dipole interaction both for 
a commensurate and an incommensurate filling of the 
lattice potential, around a small density n = 1/4. In 
the commensurate case, we observe a superfluid to in¬ 
sulating quantum phase transition. Most notably, in¬ 
troducing defects in the presence of the weak periodic 
potential allows for the realisation of defect-induced su¬ 
persolidity as originally proposed by Andreev and Lif- 
shitz [10]. In our case, the triangular crystalline struc¬ 
ture for the dipoles is either due to direct strong dipo¬ 
lar interactions or to ” pinning” of the weakly interacting 
superfluid by the lattice potential (without opening of 
a spectral gap), depending on the system parameters. 
We present a zero-temperature phase diagram which fur- 
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Figure 1: (a) Sketch of setup: Particles are prepared in a two- 
dimensional self assembled crystal with an additional quasi- 
commensurate optical lattice with varying strength A. Par¬ 
ticles are colored by their numbers of nearest neighbors AV 
Interstitials : (bl) Classical steady state of two interstitials for 
A = 0. In the absence of an optical lattice, defects collapse 
to dislocation pairs which are attractive (arrows) and form a 
stable pair in a defect string [29 . Dots indicate the position 
of the undistorted lattice; (b2) For a deep lattice A = 6 with 
spacing ao the interaction between interstitials is purely re¬ 
pulsive; (b3) Interstitials for lattice spacing ao/2 and A = 6 
can form a bound pair with triangular order. Vacancies : (cl) 
Vacancies in a self-assembled crystal and A = 0 form defect 
strings, similar to interstitials [29]; (c2) Vacancies in a deep 
lattice with A = 6 and spacing ao (c2) and ao/2 (c3) are 
purely repulsive. 


nishes a complete description of the superfluid-supersolid 
and crystal-supersolid phase transition. We find that 
superfluity and gimsi-condensate fraction are compar¬ 
atively robust against finite temperature, and estimate 
possible experimental parameters for experiments with 
ground state polar molecules [321439] and weakly bound 
molecules of strongly dipolar magnetic atoms (401143] . 
These results differ qualitatively from their deep-lattice 
counterparts, where a superfluid state was found for these 
low densities [30l l3T] . 

We note that investigations of the quantum mechani¬ 
cal phases of lattice Hamiltonians in the continuum limit 
have so far focused on contact interactions [44H461. This 
has led to the prediction of generalized superfluid-Mott 
insulator transitions in continuous space 03 . as well as 
of solid and superfluid phases of, e.g, He, He 2 and D 2 
adsorbed on solid-state surfaces such as, e.g., graphite 
[45], graphane [49H52j , Alkali substrates [55] 54 . The 
present work takes a step in a similar direction for the 
case of finite range interactions. 

The methodology presented in this paper introduces a 
new toolbox for the manipulation of complex matter, in 
analogy to the techniques developed to tune and shape 


the direct interactions between particles in systems as 
diverse as classical colloids [551162] , as well as atomic and 
molecular systems in the quantum regime [44, 163H67] , 
which is the basis for the success in the realization 
of many-body phases in these systems [68M70] . The 
proposal is based on self-assembled 2D crystals of polar 
molecules. For details on the required parameters for 
trapping and interaction strengths see Ref. 03 - 

The paper is organised as follow: in Section[H]we intro¬ 
duce the model Hamiltonian. In Section III we show and 
discuss in detail results concerning defect interactions in 
the classical regime. The study of the quantum phases is 
proposed in Section IV. Section V discusses the influence 
of temperature on supersolidity as well as parameters for 
possible physical realisations of self-assembled crystals in 
cold dipolar quantum gases (44]. Finally, in Section VI 
we draw some conclusions and future outlooks. 


II. MODEL 


We consider ultracold polar molecules which are 
trapped in two dimensions with an additional triangular 
optical lattice [Fig. [lja)] described by the Hamiltonian 

H = ^ E p? + E v » + x E y)- (!) 

i i<j i 


Here m is the mass, while represents the single particle 
momentum. In the presence of a perpendicular electric 
field, the dipole-dipole interaction Vij = To/r 3 is purely 
repulsive [44 with strength To- The last term of Eq. 0 
represents the external potential of a triangular optical 
lattice with depth A 


U(x,y) = — sin 7r(x + VSy)/ao 


1 2 


( 2 ) 


n 2 


7r (—x T \Z~3y^J /uq — sin [27rx/ao] 2 , 


which can be implemented with two standing laser beams 
m- In the following energies and distances are given 
in units of T = To/ag, ao being the optical lattice’s 
constant. The particle density p is given in units of 
N/(Aoq\/ 3/2), where A is the size of the system and N 
the particle number. As a consequence, the filling frac¬ 
tion is n = 1 for particles with average distance a = ao 
and n = 1/4 for a = 2ao, respectively. We consider a 
crystal commensurate with respect to the lattice if the ra¬ 
tio between density and filling fraction is an integer value 
and the crystal has triangular symmetry. Therefore, both 
n = 1 and n = 1/4 are commensurate configurations. 

Hamiltonian Eq. 0 can be realized with cold po¬ 
lar molecules (e.g. KRb [33] . RbCs [72] . NaK [38^ or 
LiCs [37]). where particle densities are usually p <C 1 [73. 
Alternatively, Rydberg-dressed atoms (63] E7] 1741178] or 
ground state magnetic atoms [40H43] may be used, with 
p ~ 1. In the latter case, recent experiments with weakly 
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Figure 2: Effective interaction V e ft between interstitials and 
vacancies in a dipolar crystal as a function of the distance r for 
various A and r = 15 from a classical Monte Carlo simulation, 
(a) For A = 0, V e s for interstitials is purely attractive (black). 
With increasing A and lattice filling n — 1, the interaction 
can be tuned from attractive to purely repulsive. The dashed 
line is the analytical prediction A F/ksT — Y/r 3 for a lattice 
model, (b) V e s between interstitials can also be tuned with A 
in a lattice with filling n — 1/4 (which corresponds to a lattice 
spacing ao/2). For A < 6 the shape of the interaction is similar 
to the filling n = 1. For A ~ 6, however, a bound defect 
configuration becomes meta-stable [compare with Fig. [ljb3)] . 
(c) Tunability of vacancy interactions for lattice filling n — 1. 
V e ff changes sign with increasing A. (d) Tunability of vacancy 
interactions for lattice filling close to n = 1/4. Note that 
for vacancies and A > 6 hopping is basically suppressed. For 
interstitials, due to the pair interaction a sampling up to A = 
10 is possible. Note, that the roughness of the interaction is 
an artefact of the projection of a two-dimensional interaction 
to ID. The two-dimensional potential surface for A = 0 is 
shown in Ref. [29] . 


bound Er 2 molecules composed of two highly-magnetic 
ground state atoms open the way to the realisation of 
cold quantum gases with comparatively large dipole mo¬ 
ments and high densities p ~ 1 dEo]. 


III. DEFECT INTERACTIONS 

Interstitials (vacancies) are point defects which result 
from adding (removing) a single particle to (from) a self- 
assembled crystal. Due to entropy, a finite number of de¬ 
fects exists even in thermal equilibrium [3]. These point 
defects induce a long-range displacement field mm with 
a non-linear short range part that is responsible for the 
complex dynamics of isolated defects [81] and the inter¬ 
action between them [5, 29] . On a quantitative level, 
effective interactions can be understood as the change 
in free energy, AF, as a function of the defect distance 


|r| = r for various A and fixed parameters (T, n, a o) with 

E eff (r) = A F(r)/(k B T) = - ln(5(r[x] - r)). (3) 

Here, (•) denotes ensemble averages and P(r) = (£(r[x] — 
r)) is the probability to find two defects at distance r. 
The first equality corresponds to the so-called reversible 
work theorem [82] . 

The parameter r in Eq. © describes the distance be¬ 
tween the two defects and is determined using the fol¬ 
lowing protocol, already introduced e.g. in Refs. [8, 29]. 
In each time-step, a virtual triangular lattice with lattice 
spacing equivalent to the average particle distance is con¬ 
sidered. Then, each particle is associated with the closest 
virtual lattice site. This implies that in the presence of 
an interstitial there will be exactly one doubly-occupied 
virtual lattice site. The position of this lattice site is de¬ 
fined as the position of the interstitial. For vacancies, 
the unoccupied lattice site is the position of the vacancy. 
The free energy profile of two defects as a function of 
the distance r corresponds to the effective interaction 
between the two defects. In the following, Eq. © is eval¬ 
uated using Monte Carlo sampling in combination with 
the self-consistent histogram method (see Ref. [83]). The 
effective force between the defects is then the negative 
slope of the free energy. 

In free space (A = 0) vacancies and interstitials at¬ 
tract each other in all combinations (vacancy-vacancy, 
interstitial-interstitial) and form string-like defect clus¬ 
ters, as shown in Ref. [8]. Examples of this behaviour 
are given in Figs. [T] (bl) and (cl) for two interstitials and 
two vacancies, respectively. The resulting effective in¬ 
teraction potential V e s is purely attractive and increases 
monotonically with r for all cases, as shown in Fig. §a-d) 
for A = 0 (black continuous line). 

We find that the presence of an additional optical lat¬ 
tice [Fig. [lja)] changes dramatically the energetics, dy¬ 
namics and interaction of defects. In particular, by in¬ 
creasing the lattice depth Wff can become repulsive. Ex¬ 
ample results for V e s(r) between interstitials and vacan¬ 
cies for various lattice depths are shown in the upper and 
lower panels of Fig. [2] respectively, for two choices of par¬ 
ticle densities n = 1 [panels (a) and (c)] and n = 1/4 
[panels (b) and (d)]. In all cases, the figure shows that 
for interstitials the turning point where Wff turns first 
from attractive to repulsive is A ~ 1, while for vacancies 
a larger depth of A ~ 4 is required. For comparatively 
large lattice depths (e.g., A > 6) the interaction starts to 
approach the black dashed lines, which correspond to an¬ 
alytical results from a discrete lattice model introduced 
below. The sign and strength of the effective interac¬ 
tions is however density dependent for intermediate lat¬ 
tice strengths. For example, in the case of interstitials 
V/ff displays a non-monotonic dependence on A for den¬ 
sities close to n = 1/4 [panel(b)]. For the dynamics of 
defects, this implies phase separation and cluster forma¬ 
tion for 2 < A < 6. An example of this behaviour is 
given in Fig. [l] (b3). Finally, we find that for A oc 
the dynamics of defects is effectively frozen, as the en- 
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ergy necessary to hop from one site to the next becomes 
increasingly prohibitive. 

The limit of a large lattice depth in the commensurate 
crystal (with any filling) can be understood from the fol¬ 
lowing simple lattice model. In this model, particles are 
fixed to lattice positions and can hop between sites with 
given (temperature-dependent) rates. Each lattice site 
can be un-occupied, occupied or doubly-occupied. As 
above, we assume that the direct interaction between in¬ 
dividual particles is Vij = T/r 3 . We remove the diver¬ 
gence at r = 0 by fixing the energy of a doubly-occupied 
site to U. 

Let us first consider the interaction between two in¬ 
terstitials in this model. In this case, all sites will be 
occupied and two sites are doubly-occupied, for the case 
of unit filling n = 1 plus two additional particles. The en¬ 
ergy of the system is then En = 1/2 +2 Vij+2U. 

The first term corresponds to the sum of all interactions 
in the system with TV + 2 particles and the second term 
is the 2 U offset from the two interstitials. The first term 
is a function of the distance between the additional par¬ 
ticles, say TV + 1 and TV + 2. The effective potential 
between the defects reads V G s(r) = E\i(r —* oo) —Eu(r). 
The 2 U term and all contributions in the first term up to 
last term Vn+i,n+2 cancel. This is identical to the inter¬ 
action between two particles and therefore, the effective 
potential is V/v+i,at +2 = Ves = T/r 3 . 

For vacancies, the situation is less obvious. In this 
case, the system contains TV — 2 occupied lattice sites 
and two un-occupied lattice sites. Here, the distance 
between the two vacancies is defined as the distance 
between the two un-occupied lattice sites. Considering 
again V g q (r) = Eyy(r — > oo) — Eyy(r), we find that 
when summing up all energy contributions in the system 
the energy Eyy(r) behaves exactly as V e ^(r) = V^-(r). 
Therefore, also two vacancies in the lattice are repulsive 
with T/r 3 . These exact results are depicted in Fig. [2] as 
dashed lines. 

Based on this microscopic classical model of tunable 
defect interactions, in the next Section we present results 
on defect-induced quantum phases using exact quantum 
Monte Carlo simulations for bosonic dipoles. 


IV. QUANTUM PHASES 



X 


Figure 3: (Colour online) Superfluid fraction (top panel) and 
Structure factor (lower panel) as a function of the lattice 
depth A with ra ~ 7.5 for filling n=l/4, i.e. without de¬ 
fects (Nd — 0, circle), and for n > 1/4, i.e. with defects: 
Nd — 4 (full square, interstitial density 0.02), and Nd = 8 
(open square, interstitial density 0.04). 


higher densities as n > 1/3 [T71 13H I84U86] . This regime 
of lower densities is thus essentially unexplored from the 
point of view of the investigation of quantum many-body 
phases. In the following we will be specifically concerned 
with the emergence of supersolid behaviour for densities 
close to the commensurate filling n = 1/4. 


In our analysis, we use an exact numerical quantum 
Monte-Carlo algorithm in the continuous-space path in¬ 
tegral (PIMC) representation [44j j87j US]- Our PIMC 
code is based on the so-called worm algorithm, which is 
known to efficiently provide numerically exact estimates 
of thermodynamic quantities such as the superfluid den¬ 
sity and the structure factor, which can be used as or¬ 
der parameters to determine the nature of superfluid and 
solid phase, respectively. The superfluid density reads 


fa = 


rnksT 

Wn 


<w 2 >, 


(4) 


with w = ( w x ,w y ) the winding number estimator along 
the x and y directions [89]. The static structure factor is 
instead defined as 


The present section examines the applicability of 
classical predictions to the quantum regime, considering 
commensurate and incommensurate filling values around 
n = 1/4. We chose this value as the classical groundstate 
configuration at n = 1 /4 is the commensurate triangular 
crystal that best minimises lattice-induced frustration: 
the free-space triangular lattice is very little distorted 
by the external lattice potential. In addition, previous 
studies using tight-binding models for bosonic particles 
valid in the limit of very large lattice depths have mainly 
focused on experimentally challenging (for molecules) 


$(*) = ^CE e ~ ik{ri ~N> ( 5 ) 

ij 

with |k| = k a crystal vector, and characterizes diagonal 
order. We use up to TV=188 particles and about 750 sites 
to minimize finite-size effects and defect concentrations 
of up to four percent. 

Quantum phases of Eq. 0 with A = 0 have been 
investigated for the case of bosonic polar molecules in 
Refs. [44b46l [90] . The phase diagram is characterised 
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Figure 4: (Colour online) Phase diagram rd vs. A for Eq. (1). 
f 8 and S max for rd = 15 (dashed line) are reported in Fig[5j 
see text. Error bar on rd has been estimated using a fine size 
scaling analysis. 


by a melting quantum phase transition from a triangu¬ 
lar crystal phase [with S max ^ 0 and f s = 0] to a ho¬ 
mogeneous superfluid [with S rnax = 0 and f s ^ 0] by 
decreasing the interaction strength, defined as 


Td 


Tom 
ah 2 ’ 


(6) 


of the so-called pinning quantum phase transition that 
has been predicted in one dimension [79] E2 |93]. For 
short-range interactions, the one-dimensional analog has 
been observed experimentally in ID [92l EU while the 
theree-dimensional case has been studied via quantum 
Monte-Carlo simulations in Ref. 03 - 

The scenario of the superfluid-crystal quantum phase 
transition described above changes drastically when de¬ 
fects are introduced: in this case f s remains finite but 
not unitary, as expected for a non-homogeneous super¬ 
fluid [95] . This is shown in Fig. [3] (upper panel)] in the 
parameter regime 8 < A < 12, where the superfluid den¬ 
sity increases with the concentration of interstitial defects 
(see, in particular, the cases A = 10 and 12). 

Interestingly, we find that in the whole parameter 
range A > 8 the order parameter S max is finite and 
maximal at the crystal wave-vector k = (47r/3ao,0) cor¬ 
responding to the commensurate dipolar crystal with 
n = 1/4 (as in the case described above), implying a 
diagonal crystalline order with a periodicity that is dif¬ 
ferent from that one of the underlying lattice potential 
[96] . This is in contrast, e.g., to the case with short-range 
interactions m, where the periodicity of the (Mott) in¬ 
sulating phase is trivial, in that it coincides with the one 
of the underlying lattice. 

Here, the coexistence of a finite S max and f s > 0 
in the parameter regime with A > 8 demonstrates the 
realisation of a supersolid state of matter with coexisting 
diagonal order and superfluidity. 


below a critical value r c d = (18 db 3). Here we are inter¬ 
ested in investigating the phase diagram as a function of 
rd and of the depth A of an additional optical lattice. 

Figure |3] (upper panel) displays the ground state limit 
of f s vs. A for an interaction strength rd = 7.5 < r^, 
corresponding to a superfluid for A = 0. In the figure, 
we consider both the case of a lattice commensurate with 
the dipole density at filling n = 1/4 (black circles) and 
the case of non-commensurate filling with a small density 
of interstitial defects (equaling 0.02 and 0.04 for the red 
empty and full squares, respectively). 

For n = 1/4 (black circles), the figure displays a sud¬ 
den drop of / s , suggesting a quantum phase transition 
from a homogeneous superfluid (f s = 1) to an insulating 
phase (f s = 0) induced by the increase of the strength 
of the lattice potential, occurring at the critical value 
A « 8. This picture is corroborated by a sudden increase 
of the static structure factor at the same value A ~ 8 
[Fig- H (lower panel)], showing that the insulating phase 
is in fact a crystal. This transition is thus driven by the 
suppression of quantum kinetic energy with increasing 
A [91 . As we discuss below, for physical realizations 
with, e.g., polar molecules KRb or Er 2 , the value of A ~ 8 
corresponds to a very small value of the lattice depth, 
where the two-dimensional band structure is not formed. 
Thus, this quantum phase transition can be regarded to 
happen in the continuum, as a two-dimensional analog 


Figure [4] shows the complete zero-temperature phase 
diagram as a function of and A, again keeping the de¬ 
fect density fixed at 0.02. As described above, we find 
regions of superfluid, crystal and supersolid behaviour. 
In particular, for A = 0 we re-obtain the phase diagram 
for a bosonic dipolar gas discussed above [33], charac- 



Figure 5: f s and S max vs. A for N = 26 (circle) and N = 52 
(square) considering rd = 15 (see dashed line in Fig. [4|. 
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terised by the quantum melting transition of the trian¬ 
gular crystal into a homogeneous superfluid at r d ~ 18. 
We find that increasing A from 0 initially has the simple 
effect to shift the quantum melting transition to smaller 
values of < r d , consistent with the example of Fig. [3] 
(black empty dots). However, the nature of the transi¬ 
tion changes dramatically above A > 5: a defect-induced 
supersolid phase intervenes between the superfluid and 
the crystal. For A > 6 all phases can be observed by 
simply tuning r^. 

An interesting example of this latter situation is shown 
in Fig. [5j which presents results for = 15, correspond¬ 
ing to a cut in the phase diagram of Fig. [4] along the indi¬ 
cated dashed line. Figure [5] shows that, while S' max grows 
monotonically with increasing A > 2, the superfluid den¬ 
sity f s displays an interesting re-entrant behavior, it is 
f s 1 for A < 2, then drops to zero for A 4, and then 
becomes again finite with f s ~ 0.2 - 0.4 for A > 6; this 
latter behavior corresponds to a supersolid. We notice 
that in this parameter region with < r c d the crystalline 
structure is purely imposed by the presence of the lattice 
potential, which pins the strongly interacting superfluid, 
however, crucially without opening a gap. This is similar 
to the lattice supersolid discussed in Refs. on no, how¬ 
ever it occurs for shallow lattice depths, where the band 
structure is not formed (see below). 

For rd > r d the triangular crystalline structure is 
present also for A = 0. As the crystal is essentially clas¬ 
sical (however, see [39]), we expect that the results on 
the tunability of defect interactions derived in the Sec.Ill 
above should provide directly insights into defect dynam¬ 
ics in this parameter regime. Indeed, in the PIMC quan¬ 
tum calculations we find that the effect of a sufficiently 
deep lattice in this crystalline case (e.g., A > 8) is to origi¬ 
nate a finite superfluid density coexisting with crystalline 
order, when a finite density of defects is present. This 
defect-induced supersolidity is only possible for effective 
repulsive interactions between the defects, as would be 
predicted by the classical results given above. 

Finally, for sufficiently large A > 11 the superfluid 
fraction vanishes altogether, and the groundstate evolves 
into an insulating lattice-type crystal. This is similar to 
the observed frozen dynamics in the classical regime. 

In order to visualise the difference between the solid 
and supersolid phases, Fig.[6]shows snapshots of the pro¬ 
jection of world-lines onto the xy-plane taken from the 
PICM simulations, obtained by tracing over imaginary 
time, at n = 1/4 in the absence of defects [Fig. |6][a)] 
and with two interstitials [Fig. [6](b)] for = 15 and 
A = 10. As explained in literature [88], these projections 
(that are for illustration purposes only) are the closest 
representation of the square of the wave-function for the 
many-body system that can be obtained in a simulation, 
where overlapping paths imply exchanges among the 
bosonic particles and superfluidity. The figure shows 
that in the absence of defects, paths remain localised 
around the local minima of the lattice potential U(x,y). 



Figure 6: (Colour online) Quantum Monte-Carlo snapshots 
with iV=80 (a) and N= 88 (b) particles with rd ~ 7.5 (320 
sites, square points) and A=10. In the presence of interstitials 
(right panel) delocalized paths emerge (red thick line) and 
result in a finite f s . (c) Spherical averaged one-body density 
matrix (with same set of parameter of panels a and b) for 
n— 1/4 (crystal, open symbol) and for a supersolid phase on 
lattice (full symbol). The error bars lie within point size. 

However, in the supersolid phase localised paths coexist 
with paths that are delocalized throughout the system, 
representing cyclic exchanges (permutations) among 
bosons [see Fig. §b), red thick line]. This residual 
exchange mechanism is consistent with defect-induced 
supersolidity as originally proposed by Andreev and 
Lifshitz HO] and only recently demonstrated via exact 
theoretical techniques for bosons with cluster-forming 
interactions m ■ Our analysis shows that defect-induced 
supersolidity can be originated in the continuum also for 
non-cluster forming liquids, using periodic external po¬ 
tentials. Apart of the qualitative difference, as explained 
below, this should provide quantitative advantages in 
the experimental realisation of the supersolid phase, as 
it could results in, e.g, higher temperatures than possible 
in the tight-binding regime. 

Particle delocalization is also reflected in the (quasi-) 
condensate fraction, which is easily accessible in experi¬ 
ments m, defined as the asymptotic (i.e., r —> oc) be¬ 
haviour of the angle averaged one-body density matrix 

n 0) = ^ J dfl j drn(r,r'), (7) 

with 


n(r,r') = (^(r )^ f (r')), 


(8) 
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and -0(r) [^(r)] the particle annihilation [creation] op¬ 
erators at position r m • In the presence of long-range 
off-diagonal order associated to a finite condensate frac¬ 
tion, n(r, r') factorizes at large separation |r — r'| as 

)V' t ( r ')) -> <K r )<K r ')> (9) 

with (j){ r) the condensate wave-function. Employing the 
same set of parameters of Fig. [6^b) , Fig. [HJc) shows that 
a constant value for the supersolid phase with finite 
defect concentration is here realised (full squares), corre¬ 
sponding to a finite (quasi -)condensate fraction at T ~ 0. 
The latter disappears in the case of commensurate filling 
n = 1/4 [parameters as in Fig. [6](a)] , where n(r, r') 
decays exponentially with distance [empty square in 
panel (c)], as expected for an insulating crystalline phase. 


V. FINITE TEMPERATURE AND POSSIBLE 
PHYSICAL REALIZATIONS 

The supersolid phase described here is considerably 
resilient towards finite temperature effects. An example 
of this is shown in Fig. [7| where we present results 
for both f s and S rnax as a function of the rescaled 
temperature t = ksT / (h 2 n/m), expressed in units of 
the quantum kinetic energy at the mean inter particle 
distance h 2 n/m. The figure shows that supersolid 
behaviour survives up to temperatures of the order of 
t = ksT / (h 2 n/m) <1. The latter is consistent with 
a BKT transition for f s with a comparatively large 
transition temperature Tbkt (see below). 



t 


Figure 7: (Colour online) Superfluid fraction (top panel) and 
Structure factor (lower panel) of the reduced temperature t m 
ksT / (h 2 n/m) with A=10 and Nd — 4 (full square, interstitial 
density 0.02), and Nd = 8 (open square, interstitial density 
0.04). 


A. Polar molecules: KRb 

As an example of experimental realisation of the phase 
above, here we consider first a gas of bosonic 39 K 85 Rb 
molecules with dipole moment d ~ 0.5 Debye trapped 
on a lattice with spacing ao = 400 nm. The lattice 
recoil energy is Er/H = h/fima^) ~ 2.49 kHz (in fre¬ 
quency units), while the unit of energy T = Yo/a^ reads 
Y/h ~ 0.59 kHz. By re-expressing the value A = 10T in 
convenient Er units, the numbers above suggest that the 
supersolid phase described above would be realised for a 
weak lattice potential with a depth of just > 2.3 Er. This 
corresponds to a situation where the band structure of 
the two-dimensional lattice has not formed, which is con¬ 
sistent with our claim that this quantum phase transition 
occurs in the continuum. The quantum kinetic energy at 
the average particle density E^ n /h = Tm/(27rraKRb) i n 
the absence of the optical lattice reads E^[ n /h ~ 130 Hz. 

. Thus, from the results of Fig. [7] we obtain that the su¬ 
persolid phase would survive up to temperatures of the 
order of T c ~ 6 nK or larger. 


B. Magnetic quantum gases: Er2 

Magnetic quantum gases of Er have been recently 
trapped in lattices with a small spacing ao = 266nm [98] . 
For Er 2 molecules with magnetic dipole moment of 
Fr 2 = 14 i±b (with /i b the Bohr magneton) this implies 
Y/h ~ 0.14kHz, which is comparable to the case of polar 
molecules. Extrapolation from our numerical results (see 
Fig.g implies that a lattice depth of order of A ^ 10T 
is necessary to induce supersolidity. 

Here, the recoil energy is Er ~ 2.1kHz, and thus the 
condition on the lattice depth A for inducing a super¬ 
solid behaviour reads A ~ 10T ~ 0.7 Er, which is well 
in the continuum limit. The quantum kinetic energy 
reads E\^ n /h = hn/(27rmEr 2 ) — 105Hz, for a low den¬ 
sity n = l/(4ag). This corresponds to a critical tempera¬ 
ture for observing supersolidity of the order of T c > 5nK, 
which is within experimental reach. We note that for the 
case of strongly dipolar magnetic gases, reaching den¬ 
sities of order unity in the lattice is possible. This in 
principle would allow for the observation of some of the 
effects described above, such as the superfluid-insulator 
quantum phase transition for commensurate lattice fill¬ 
ings, at considerably higher temperatures. 


VI. CONCLUSIONS AND OUTLOOK 

The results presented here may help to settle the 
long-standing question of the role of defect interactions 
for the formation of the supersolid phase da cm] and 
enable the experimental realization of a supersolid with 
quantum dipolar gases in combination with tunable 
optical lattices. We note that while here we have 
focused on specific examples with fixed defect densities, 








defect-induced supersolidity is expected to appear for a 
large range of densities for a careful choice of interaction 
parameters, e.g., similar to the case of Ref. [115]. This 
should make the observation of defect-induced super¬ 
solidity possible in experiments where particle density 
is tunable within a few percent only. The tunability of 
defect interactions opens also interesting prospects for 
the observation of other phases, such as solitons and 


breathers in the classical and quantum regimes [79] [99]. 
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